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We investigate the mass differences of charm and bottom heavy baryons, using MILC lattice 
gauge configurations with 2+1 dynamical sea quark flavors. We extend our previous work to three 
lattice spacings: fine (a i=s 0.09), coarse (a « 0.12), and medium-coarse (a « 0. 15) ensembles. For 
extrapolations and interpolations, we apply simultaneous quadratic fits and simple linear fits with 
the full QCD data points for the valence quark masses as well as the sea quark masses. 
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1. Introduction 

Heavy baryons have been investigated by both experimental and theoretical approaches. From 
experiment, the singly charmed heavy baryon mass spectrum is well known; however, the other 
heavy baryon masses are only crudely known. Recently, mass measurements of the singly bottom 
baryon have been published by the D0 [|l|] and CDF [Q] collaborations. From lattice QCD, there 
are several quenched calculations [3-8] and dynamical sea quark flavor simulations [9-1 1] for the 
heavy baryon mass spectrum, and most results are in fair agreement with observed values. In this 
work, we extend our previous dynamical sea quark flavor simulation []9|][10] with a larger data set 
and more systematic analysis methods. 

We use the MILC fine (a m 0.09), coarse (a 0.12), and medium-coarse (a w 0.15) lat- 
tices [0|. We apply two local interpolating operators []3|] and construct heavy baryon two point 



functions with improved clover heavy quark propagators with the Fermilab interpretation []13j] and 
improved staggered light quark propagators. We extract the mass differences from the ratio of 
propagators and extrapolate to the chiral limit using simultaneous quadratic fits and simple linear 
fits with the full QCD data points. 

2. Formalism 

We used two types of interpolating operators ^ for the singly heavy baryons, which are 

5 = E ahc (y a l T Cy 5 y^ H and ^ = e ahc {¥? Cy^)^ , (2.1) 

where e a b c is the Levi-Civita tensor, \\f\ and \j/2 are light valence quark fields for up, down, or 
strange quarks, is the heavy valence quark field for the charm or the bottom quark, C is the 
charge conjugation matrix, and a, b, and c are color indices. Basically, is the operator for 
s n = + and G» is for s n = 1 + , where s % is the spin parity state of the light quark pair. Therefore, 
<^5 describes total spin J = \ baryons (A#, E#), and ffp describes J = A baryons (£#, E' H , £2#) as 
well as J = | baryons (L H , E' H , Cl H ). For doubly heavy baryons, we can simply interchange the 
heavy quark field and the light quark fields in Eq. O] [Bj. 



We apply the method of Wingate et al. [14] to combine staggered propagators for the light 



valence quarks and a Wilson type (clover) propagator for the heavy valence quark. Since we use 
staggered light quarks, we needed to consider taste mixing. In our previous work [|K^[|l5|], we 
found that there is no taste mixing between the light quarks. However, due to cancellations among 
the copy indices, we cannot separate spin J = i and J = 3 states from the operator using 
standard spin projection operators. 

3. Data analysis 

Since we use the Fermilab inteipretation Jl3| ] for the heavy quarks, the absolute mass from the 
simulation is not a physical quantity. The physical mass M p i, y of the baryon is 

M phy =M ca , + A, (3.1) 
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where M ca i is the simulation result and A is a constant mass shift. We can calculate the constant 
mass shift from calculations of the kinetic mass with non-zero momenta or heavy-light meson 
spectroscopy. In this project, however, we estimate mass differences, which are directly measurable 
quantities without determining A. 

Since we are interested in mass differences, we can take the ratio of two propagators and fit 
the mass difference directly. A fit model function P r (t) of the ratio of the propagators is 
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(3.2) 



Since the propagators are calculated using the same gauge configurations, the propagators in nu- 
merator and denominator are correlated with each other. Using this method, we expect smaller 
statistical errors and more stable fits. 

We can take the ratio of the propagators, for example, between E' c and A c , 



Crit) 



(3.3) 



where k is the configuration index. The covariance matrix of the ratio is estimated using bootstrap 
sampling. 

We display an example of this method for E c and A c on the fine lattice with dynamical quark 
masses mi = 0Am s and 0.0062 light valence quark mass in Fig. |[ The jittery pattern appears 
strongly in the ratio (Fig. [I] (a)), because the opposite parity state contributions are substantially 
different for numerator and denominator. In fact, E' c contains large opposite parity state contribu- 
tions, while Ac contains very weak opposite parity state contributions. 

We obtained the most reasonable fit from a fit model function with one non-alternating phase 
exponential and one alternating phase exponential. That fit is shown in Fig. [j](b). The maximum 
distance of the fit is D max = 22 and a plateau appears where D m ; n > 8. 



-a o.i 

3 



I ' ' ' ' I ' ' ' ' I ' ' 
(a) H' - A c / 0.4m^ fine, m0.0062 



□ □ 



mm 



mm 



<D 0.2 



to 
S 



1 ' ' ' ' I ' ' 
(b) E' c - A„ / 0.4m fine, m0.0062 



□ 
□ 



10% CL 
50% CL 
90% CL 



[] 



'[Hpip[]]5 C] n 



[][] 



10 20 

time 



Figure 1: The ratio of the propagators (a) and the two particle fit of the ratio (b). The size of the symbols on 
(b) is proportional to the confidence level (CL). 
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The mass difference can also be obtained by first fitting the propagators individually, and then 
subtracting the masses to obtain the difference. We can compare this method shown in Fig. || to 
the ratio method. As we expected, the individual mass fits for E' c and A c (Fig. ||) are worse than 
the fit of the ratio (Fig. [I] (b)). For E' c , the confidence levels are fine; however, the fluctuation 
of the plateau is larger (the systematic error is larger), and for A c , fits are not stable. From this 
comparison, we conclude that the ratio method is better. 
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Figure 2: The individual fits for E'. (a) and A £ - from the propagators. The square symbols indicate successful 
fits and the cross symbols indicate unstable fits. The size of the symbols is proportional to the confidence 
level. We obtained these fits from a fit model function with one non- alternating phase exponential and one 
alternating phase exponential. 

Once we obtain the mass differences, we need to perform extrapolations to the physical quark 
mass. In this analysis, we extrapolate the valence quark masses and the sea quark masses simulta- 
neously, using a quadratic fit model function given by 

^quad = Co + c x m\ + c 2 mj + c 3 m s + c 4 m sea (3.4) 

where cq to C4 are the fitting parameters, m/ is the light valence quark mass, m s is the strange valence 
quark mass, and m sea is the light sea quark mass. We fixed C3 = except for the coarse lattice, since 
only the coarse lattice has multiple (three) strange valence quark masses. We would prefer to use 
a heavy baryon mass formula based on partially quenched chiral perturbation theory, but do not 
know of one. We can also perform the extrapolations with only the full QCD data points, i.e., the 
points for which the valence quark masses match the corresponding sea quark masses. In general, 
the errors of the full QCD extrapolations are larger than those of the simultaneous quadratic fits, 
because the simultaneous fit uses more data to constrain the fit parameters. 



4. Results and discussion 

Parameters of the ensembles [12] are summarized in Table [jj. The scale of each ensemble was 
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a (fin) 


ami/am s 


L(fm) 


P 


size 


n/a 


# confs 


~0.15 


0.0097 / 0.0484 


2.4 


6.572 


16 3 x 48 


2.1356 


631 


~0.15 


0.0194 / 0.0484 


2.4 


6.586 


16 3 x 48 


2.1284 


631 


~0.15 


0.0290 / 0.0484 


2.4 


6.600 


16 3 x 48 


2.1245 


440 


~0.12 


0.007 / 0.05 


2.4 


6.76 


20 3 x 64 


2.6349 


545 


~0.12 


0.01 / 0.05 


2.4 


6.76 


20 3 x 64 


2.6183 


590 


~0.12 


0.02 / 0.05 


2.4 


6.79 


20 3 x 64 


2.6477 


452 


~0.09 


0.0062 / 0.031 


2.4 


7.09 


28 3 x 96 


3.7016 


534 


~0.09 


0.0124/0.031 


2.4 


7.11 


28 3 x 96 


3.7138 


557 



Table 1: MILC lattice parameters. The three nominal lattice spacings in the first column (0.15, 0.12, and 
0.09 fm) correspond to medium-coarse, coarse, and fine ensembles, respectively. The bare light (strange) 
sea quark mass in lattice units is ami (arris)- The spatial size of the lattice is L and j3 = IQ/g 2 , where g is the 
bare gauge coupling. In the sixth column we show r\/a calculated from a global fit to r\ values on all our 
ensembles [jig], 

determined by a length scale r\ from the static quark potential using f n as an input parameter JT7| ] : 

n =0.3108(15) (i^)fin. (4.1) 



We present mass splittings in MeV for singly heavy baryons in Fig. ||. We investigate five 
J p = 2 + singly heavy baryons (Ah, E#, E' H , £Ih)- There are ten possible mass splitting 
combinations. Since not all of the combinations are independent and all dependent combinations 
are consistent each other, we display only four independent combinations on the figure. We display 
the comparison of the singly charmed baryons and the PDG data [ |l8| ] in Fig. || (a). The error bars 
of our calculations are larger than those of experiment, but smaller than or comparable to those 
of previous lattice calculations [3-8]. For the singly bottom baryons in Fig. ^] (b), we compare 
with experimental results measured by the CDF and D0 collaborations and a dynamical 

calculation result by Lewis and Woloshyn [11]. Our error bars, in general, are smaller than those 
of Lewis and Woloshyn, especially for the fine lattice result. 

Since the singly charmed baryons are experimentally well known, let us examine more closely 
our singly charm results. The mass differences between S c — A c , E' c — E c , and Cl c — E' c are in good 
agreement with the PDG data [|j||; however, the mass difference £l c — A c is not. Note that A c and Sc- 
are calculated from the operator ^? 5 , and £ c , E^., and Q. c are calculated from the operator Thus, 
the disagreement occurs on the mass differences of the heavy baryons that come from the different 
operators. On the fine lattice result, especially, this pattern is quite evident. We checked the other 
mass differences with the different operator combinations, and we obtained similar discrepancies 
from the PDG data [18]. We hope that this can be resolved by studying the constant mass shift due 
to the heavy quarks or the hyper-fine structure of the singly heavy baryons. Resolving this puzzle 
is a high priority. 

We also present the doubly charmed and bottom baryons in Fig. |4j. We display our fine and 
coarse lattice results with those of Lewis et al. [Q] for the doubly charmed baryons and Lewis and 
Woloshyn [11] with dynamical sea quarks for the doubly bottom baryons. We set the scale of the 
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Figure 3: Independent mass differences of J p — | singly charmed baryons (a), and singly bottom baryons 
(b). In the legend, "full" indicates the full QCD data point fit result, and "quad" indicates the simultaneous 
quadratic fit result. 



y-axis of the figure from the other group's results, because they calculated the absolute masses, 
while we calculated the static masses. For the doubly charmed baryons, we added 193 and 438 
MeV, and for the doubly bottom baryons, we added 2.833 and 4.223 GeV to our fine and coarse 
lattice results, respectively. Our results agree fairly well with prior results. The mass differences 
between Ehh and Q.hh are about 100 MeV, which is about the strange quark mass. The hyper-fine 
splittings are about 30 ~ 80 MeV. Moreover, the hyper-fine splittings of the doubly bottom baryons 
are smaller than those of the doubly charmed baryons, because of heavy quark symmetry. 
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Figure 4: The mass spectrum of doubly charmed and bottom baryons. The error bars are statistical only. 
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